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ABSTRACT. We generalize the standard Poisson summation formula for lattices so that it operates on the level of theta series, 
allowing us to introduce noninteger dimension parameters (using the dimensionally continued Fourier transform). When 
combined with one of the proofs of the Jacobi imaginary transformation of theta functions that does not use the Poisson 
summation formula, our proof of this generalized Poisson summation formula also provides a new proof of the standard 
Poisson summation formula for dimensions greater than 2 (with appropriate hypotheses on the function being summed). In 
general, our methods work to establish the (Voronoi) summation formulae associated with functions satisfying (modular) 
transformations of the Jacobi imaginary type by means of a density argument (as opposed to the usual Mellin transform 
approach). In particular, we construct a family of generalized theta series from Jacobi theta functions from which these 
summation formulae can be obtained. This family contains several families of modular forms, but is significantly more 
general than any of them. Our result also relaxes several of the hypotheses in the standard statements of these summation 
formulae. The density result we prove for Gaussians in the Schwartz space may be of independent interest. 



1. Introduction 

Consider a lattice A C K" and a sufficiently well-behaved function F : M" — > R. [Taking F to belong to 
the Schwartz space S^(M. n ) is sufficient, and is what we shall do in our later generalization.] The standard Poisson 
summation formula then says that 

tth ^ detA P t^ 
Here A* is the lattice dual to A, dot A denotes the volume of a Voronoi cell of A, and 



F(p) := [ F(x)e- 2mx - p dx 



denotes the Fourier transform of F. (We use a tilde here so that we can reserve the circumflex for our more general, 
dimensionally continued Fourier transform.) We wish to construct a dimensionally continued version of this result. 

This problem was originally inspired by a condensed matter physics investigation involving the dimensional con- 
tinuation of electrostatic lattice sums, computed using the Ewald method (see, e.g., l20l for a modern exposition of 
this method), as described in 1211 , However, the ensuing discussion is a purely mathematical offshoot of this investiga- 
tion. For one thing, the results we were able to prove do not include the physically relevant case of a slowly decaying 
function, even though we have numerical evidence that the results still hold in this case. Nevertheless, the methods 
used here might also be applicable to the dimensional regularization of lattice sums: See (8) for an approach using 
zeta functions and the Mellin transform. 

In addition to its use in performing lattice sums, the Poisson summation formula is a prominent tool in most 
other problems involving lattices, from point counting problems in number theory (discussed by, e.g., Miller and 
Schmid ll26l ). to the mathematical theory of diffraction — see, e.g., the review by Lagarias |23| . While we do not 
consider diffraction in this article, it is possible that our results might have some application to diffraction problems, 
or that results in the theory of diffraction might provide some further avenues for generalization of the results presented 
here. In particular, the quasicrystal summation formula given in Theorem 2.9 of Lagarias Il23l seems to offer attractive 
possibilities for generalization. Additionally, Baake, Frettloh, and Grimm [2| apply an integer dimension result that is 
very similar to our formula to a problem in diffraction. 

If we specialize to the case where F is a radial function, we can obtain a dimensionally continued Poisson summa- 
tion formula in terms of the lattice's theta series — see Theorem[T]for the (particularly simple) version for 7L d (d > 1). 
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One can then generalize this result to a reasonably large family of theta series-like functions constructed out of linear 
combinations of products of powers of Jacobi theta functions; the generalization is given in Theorem[2] This family 
of generalized theta series includes all the theta series of lattices given in Chap. 4 of Conway and Sloane iflOl (except 
for the general forms of those for the root lattice Ad and its translates), as well as several families of modular forms 
(including all the classical modular forms of even weight), as discussed by Rankin BP . and more compactly stated 
as Theorem 12 in Pache |30|. However, the family we consider is considerably more general, since we only require 
integer powers to reproduce the aforementioned theta series from Conway and Sloane and the families of modular 
forms from Rankin, while here the powers can be arbitrary nonnegative real numbers. 

This relation between modular transformations and summation formulae is not new. In fact, much of the basic 
result we obtain was first given by Bochner |5 |. [The summation formula is stated in Bochner's Eq. (80) in an integral 
form and in his Eq. (129) using sums; a nice form of sufficient hypotheses on the function being summed are given in 
Bochner's Theorem 9. Also see Theorem 4.9.1 in Q for a more streamlined presentation of the summation formula 
(using sums), though with some obvious typos, and without a statement of the hypotheses concerning the function 
being summed.] Our method of proof bears some similarities to Bochner's, though Bochner uses Fourier-Laplace 
transforms as the basic analytic tool, while we use a density result. (However, Bochner's method allows for weaker 
hypotheses.) Our overall focus is also somewhat different, in that we are concerned with the summation formulae 
derived from a specific family of generalized theta series, while Bochner only considers a general modular relation. 
As far as we know, this specific construction of a family of dimensionally continued Poisson summation formulae has 
not been considered previously. While Guinand gives an explicit analogue to our Theorem[T]in Sec. 10.5 of fl31 . this 
is only stated for integer dimensions, even though Guinand cites Bochner as a reference for the formula. (Guinand 
also cites his own general summation formula from fl4*l . which is only given for integer dimensions.) Note that we 
only became aware of the work of Bochner and Guinand in the process of publication, through the reference to lfT31 in 
the review by Lagarias 11231 that the referee suggested to us. 

In more recent work, the relation has reached its most refined form in the association between automorphic forms 
and Voronoi summation formulae (see, e.g., Miller and Schmid l26l for a review of recent work). However, work on 
this correspondence first arose in the context of transformations related to the functional equation for the zeta function, 
inspired by a question from Voronoi on analogues of the Poisson summation formula — see, e.g., 11361 and references 
therein. (We call particular attention to the work of Ferrar fl"3l ; see Theorem 10.2.17 in Cohen J7] for a more modern 
discussion of a very similar result.) Additionally, Baake, Frettloh, and Grimm Q give a (distributional) radial Poisson 
summation formula in their Theorem 3 in a form that is very similar to our dimensionally continued form. However, 
they do not show how to dimensionally continue the lattice (or, indeed, mention theta functions explicitly), and their 
proof (which relies on the standard Poisson summation formula) only holds for integer dimensions. There are also 
discussions of similar formulae — these derived from modular transformations — at the beginning of Chap. 4 of Iwaniec 
and Kowalski HI 81 . and in Sec. 10.2 of Huxley IfTTI — what Huxley terms the Wilton summation formula. These 
formulae are presented in what appears to be a dimensionally continued form, though their hypotheses assume integer 
dimensions. Regardless, the summation formula in Iwaniec and Kowalski and Huxley's Wilton summation formula 
are derived from cusp forms, while our result (in the language of modular forms) does not require the vanishing of the 
constant term in the form's Fourier series (the defining characteristic of a cusp form). 

We make considerable use of various special functions in the following discussion, so we provide references for 
the conventions we use. The dimensionally continued theta series for which our summation formula holds will all 
be constructed out of Jacobi theta functions. An appropriate introduction to these functions for our purposes is given 
in Chap. 4 of Conway and Sloane fTOfl [though there is an important difference in our notation, as discussed below 
Eq. {3)]. We shall also encounter the confluent hypergeometric limit function and the Hermite functions (plus a 
brief appearance by a Bessel function). These are discussed in an appropriate way by Andrews, Askey, and Roy JT| 
(though we slightly streamline their general hypergeometric function notation, since we only consider the confluent 
hypergeometric limit function). In general, the Wolfram Functions Site BT1 is a good resource for information about 
the various special functions we employ. 

Additionally, we shall need the apparatus of tempered distributions introduced by L. Schwartz ||34| . A convenient 
overview of the properties we need is given in Sec. V.3 of Reed and Simon 1321 . (The Appendix to that section also 
provides another exposition of the iV-representation we use.) 



2. Ingredients 
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2.1. Theta series. Here we recall various facts about theta series and theta functions that we shall need for the rest of 
our discussion, following Chap. 4 of Conway and Sloane 1101 . The theta series of a lattice A is defined by 

Q A (q) :=]>>' fe ' 2 . 
fceA 

(This is often treated as a formal series, but converges for q € C, \q\ < 1.) The utility of the theta series stems from the 
fact that the coefficient of q l in the expansion of 9a (?) in powers of q gives the number of points in the intersection 
of the lattice and a sphere of radius \/l centred at the origin. Thus, if we write 

oo 

(2) e A ( g )= : ^iv,<A 

1=0 

then 

oo 

£/(!*!) =X>/(v^T). 

fceA i=o 

[Nota bene: We have written the radial function F as /(| • |), and shall only consider these radial parts in the sequel.] 
Examples of dimensionally continued theta series for families of lattices include the e?-dimensional cubic lattice Z d , 
with Q Z d(q) = tfj(q), and the root lattice D d , with <d D d(q) = [&j(q) + ■dj(q)]/2. (See Chap. 4 in Conway and 
Sloane iTTOl for further examples.) Here 

oo oo oo 

(3) d 2 {q)--=W 4 Y.i l2 ~ 1 ' *3(<z) :=i + 2£y, Mq)--=i + 2j2(-<if 

1=1 1=1 1=1 

are Jacobi theta functions (where $2 is defined for future use). 

Nota bene: It is often customary to take theta functions and theta series to be functions of a complex variable z, 
instead of the nome q = e l7TZ that we have used here. We have chosen to regard the nome as fundamental since we are 
primarily interested in the expansions of these functions in powers of q. However, when discussing transformations of 
these functions, it is considerably more convenient to regard them as functions of z. On the few occasions where we 
do this, we shall use an overbarto denote the difference, e.g., Q\(z) := @A(e I7rz ). (In the literature on summation 
formulae derived from automorphic forms, one thinks of our expansions of a in q as the Fourier coefficients of 0a-) 

Since the Poisson summation formula involves the dual lattice, we need to know how to obtain its theta series. This 
is given by the Jacobi formula [Eq. (19) in Chap. 4 of Conway and Sloane iflOl l. which states that 

(4) 0a* (z) = VdrtA(i/z) d / 2 Q A (-l/z), 

where d is the dimension of the lattice. The Jacobi formula is typically proved using the Poisson summation formula 
[see, e.g., the discussion leading up to our Eq. (8)]. However, all we need in our discussion is the intimately related 
Jacobi imaginary transformation of the Jacobi theta functions (also known as the modular identity or reciprocity 
formula for the theta functions), i.e., 

(5) M-l/z) = (zA) 1 /¥ 4 (z), h(-l/z) = {z/i) l ' 2 h{z). 

(The first of these is also true with the labels 2 and 4 switched.) The standard proof of these identities is a direct 
application of the Poisson summation formula, but there are alternative proofs that are independent of it. For instance, 
one such proof is given in Sec. 21.51 of Whittaker and Watson l40l . while Bellman's text |4'] discusses several others — 
see, in particular, Sec. 30 for Polya's derivation — in addition to the standard Poisson summation version (in Sec. 9). 
Our discussion will thus be independent of the standard Poisson summation formula (with the exception of a brief 
appeal to establish Theorem[T]for d = 1). 

2.2. The dimensionally continued Fourier transform. We also need to dimensionally continue the Fourier trans- 
form. Stein and Weiss give a dimensionally continued version of the Fourier transform for radial functions in Theo- 
rem 3.3 of Chap. IV of El, viz., 



f(p) := 2Kp-^' 2 / f{r)J [d _ 2)/2 (2vpr)r d l 2 dr 

(6) 



2n d / 2 



T(d/2) 
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(This reduces to the standard Fourier transform for a radial function when d S N.) Here the first equality gives 
the expression from Stein and Weiss (Jk is a Bessel function) and the second gives an equivalent (perhaps slightly 
neater) expression in terms of the confluent hypergeometric limit function o-Fi. The hypergeometric expression has 
the advantage of only involving one appearance of p (and being manifestly regular at p = for all d > 1), in 
addition to showing the rf-dimensional polar coordinate measure for radial functions explicitly. We shall thus use the 
hypergeometric expression exclusively in the sequel. (One can obtain the hypergeometric expression using the Stein 
and Weiss derivation — the only difference is that one uses a different special function to evaluate the final integralQ) 

For this expression to be well-defined, it is sufficient to take d > 1: One assumes d > 1 when using integral 
representations to express the result in terms of either of the two given special functions, and can also check that the 
integral is convergent for all p € M in that case, provided that / € L 1 (R+). Additionally, Eq. © reduces to the 
expected expression for d — 1 [using J-±/2(z) = \plfnz cos z or o-Fi(l/2; —z) = cos(2 v / z)]0 (Stein and Weiss 
restrict to d > 2 so that the integral they use in their derivation is well-defined, since they are only considering integer 
dimensions.) This restriction to d > 1 is necessary for other parts of our discussion, though we have numerical 
evidence that it can be relaxed. 

The following result is central to understanding why this dimensionally continued Fourier transform agrees with 
the dimensional continuation of the theta series. 

Lemma 1. For d > 1, the dimensionally continued Fourier transform (for radial functions) defined in Eq. © takes a 
Gaussian Q a {r) := e~ ar , Rea > to another Gaussian, Q a (p) = (7r/a) rf / 2 e _7r p ' a . 

Remark. Intuitively, this result follows from dimensionally continuing the well-known integer dimension result. We 
should get the same result from direct calculation using Eq. (0 since that expression was obtained using the same 
dimensional continuation procedure. 

Proof. The case d = 1 is classical. For d > 1, we use o-Fi's defining series, 

[(•) n denotes the Pochhammer symbol] and integrate term-by-term, evaluating each integral using the gamma func- 
tionQ The resulting series is the Maclaurin series for the expression we gave for Q a . The term-by-term integration 
is justified by the Lebesgue dominated convergence theorem. To see this, we use the same integral representation for 
o-Fi used in the derivation of Eq. ©, which gives, for any TV G N, 



(-7rVr 2 ) r 



< Fi(d/2;7r 2 pV) < A"cosh(27rpr) / (1 - t 2 )( d ~ 3)/2 dt, 

Jo 



where K > is a constant^ This allows us to apply the dominated convergence theorem, since the integral in the final 
term is finite for d > 1 and f Q cosh(27rpr) \e~ ar \dr is finite for Rc a > 0. □ 

Remark. The importance of this result to our discussion comes in its use in obtaining the integer dimension Jacobi 
transformation formula (and thus also the Jacobi imaginary transformations of the Jacobi theta functions) via the 
standard Poisson summation formula: For a lattice A of dimension n £ N, we have (taking Im z > so that everything 
converges) 

i / ■ \ ™/ 2 1 / • \ ™/ 2 

6 A (z) := y e-l fe l 2 = ^= ( - ) y e-^l 2 /* = -2= ( M A .(-1/*), 

which can be written as 



(8) Qa*{z) = VdetA(z/z)"/ 2 8 A (-l/z), 



'The integral representation for o-Fi we used is 07.17.07.0004.01 on the Wolfram Functions Site (41j. 
2 These identities are 03.01.03.0004.01 and 07.17.03.0037.01, respectively, on the Wolfram Functions Site (4TJ. 
3 The Maclaurin series for Fx is 07.17.02.0001.01 on the Wolfram Functions Site BTl . 
^Nota bene: We denote the set of positive integers by N, and the set of nonnegative integers by No. 
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the Jacobi transformation formula. We thus expect that the dimensionally continued dual theta series that we obtain 
using this formula will agree with the dimensionally continued Fourier transform to give a dimensionally continued 
Poisson summation formula. 

3. The dimensionally continued Poisson summation formula for Z d 
With these results in hand, we can thus write Eq. ([TJi [for a radial function F =: /(| • |)] as 



oo 



where the starred quantities come from writing the theta series of A* in the power series form given by Eq. (0, and 
we calculate / by taking the dimension parameter d to be the dimension of the lattice. (As we shall see later, what is 
important is that the d one uses here is the same d that appears in the Jacobi transformation formula.) It is clear that 
this equality holds when d £ N, by the standard Poisson summation formula. What is perhaps surprising is that the 
equality still holds for, e.g., A = Z d , with d € R (d > 1). We shall first prove the result for this simple case (Z d is 
self-dual, dct Z d = 1, and Ai = I), where it becomes 

Theorem 1. If f e J^ E (R) (i.e., f is an even Schwartz function) and d > 1, then 

oo oo 

(9) X; JVi/(v7) = 5^ 

1=0 1=0 

where the Ni are given by the power series expansion of the theta series ofL d , viz., 



(10) 9(g) = #i( q ) = 



=:£AV' 



1=0 



k=l 

and f is computed using Eq. (|6). 

However, the simplifications are primarily notational. As we shall see in the discussion in Sec. [6j the proof works 
with minimal modifications for a much larger class of 0s, including functions that cannot be the theta series of a lattice 
(even though they have an integer dimension parameter). 

Remark. The restriction that / be an even function should not be surprising: In integer dimensions, it corresponds to 
the lack of a cusp at the origin for the full radial function F = f(\ ■ |). Moreover, as Miller and Schmid note l26l . the 
standard one-dimensional Poisson summation formula is a trivial = for odd functions. 

4. A Schwartz space density result 

Since the proof proceeds by noting that the desired formula holds almost trivially for the Gaussians from LemmaQ] 
and then extends to an interesting set of functions [viz., ^ E (R)] by density, we start by establishing the requisite 
density result. 



Lemma 2. Span{x n- e Q2; "|q; > 0} = 5? (R) [i.e., the Schwartz space closure of the given family of Gaus- 
sians is all the even Schwartz functions]. 

Proof. We shall prove this by showing that 

X := Span{x i-> e - ax '\a > 0} + Span{x ^ xe~ ax2 \a > 0} 

is dense in =^(M), so its even part, Span{x n- e~ ax ~ \a > 0}, is thus dense in <5^ E (R). We shall use Corollary IV.3.14 
from Conway's text (9|, which states that a linear manifold (here X) is dense in a locally convex topological vector 
space [here J^QR)] if and only if the only element of the dual of the topological vector space that vanishes on all 
elements of the linear manifold is the zero element. 

It is most convenient to proceed by identifying J^(]R) with a sequence space, following Simon [35). (There is 
an alternative presentation of these results in the Appendix to Sec. V of Reed and Simon 11321 .) Here the sequence 
space is given by the coefficients of the Hermite function expansion of elements of =5^(R), and provides a particularly 
nice characterization of the tempered distributions [the elements of <5*"(R), the dual of J^(R)]. Namely, if a n are 
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the Hermite coefficients of / € J^(R) [i.e., a n := J m f(x)h n (x)dx, where h n is the nth Hermite function], then 
ip G J^'(R) can be written as <p(f) = J^^Lo c n a n, where c„ are the Hermite coefficients of tp, with |c„| < C(l + n) m 
for some C, m > 0. (This is Theorem 3 in Simon 13511 .) Note that Simon defines the Hermite functions to be I? 
normalized, so, we have, from the first equation in Sec. 2 of Simonfl 

h n ( X ) := ^^H n{x) , H n ( X ) := ^\ 

where the H n are the Hermite polynomials, with generating function^ 

n=0 

We can now use this generating function to show that the Hermite coefficients of x H> e~ ax2 are given by 

e ~(ax 2 +x 2 /2~2tx+t 2 ) ^ x 



*=0 



where M n := (iz 1 ' 2 2 n vX)- x f 2 is the Hermite functions' normalization factor and (3 := l/(a + 1/2). We thus have 
o-2n = A/'2, l (7r//3) 1 / 2 (/3 — l) n /nl, a2 n +i = 0, by the series expansion of the exponential. [We used Lemma 2.2 in 
Chap. VIII of Lang l24l to interchange differentiation and integration. We only need to consider the case where t 
lies in some neighbourhood of 0, so the ^derivatives of the integrand are each bounded by a polynomial in x times a 
Gaussian in x (for all t in the neighbourhood), and those functions of x are integrable over R.] Similarly, the Hermite 
coefficients of x n- xe~ ax are & 2 „ = and &2n+i = A/2 n +i(7r//3 3 ) 1 / 2 (/3 — l)™/n!. Thus, we consider 

Ep,±(x) := {P/^e-^ 2 ± (^/n^xe-^ 2 , 

which has Hermite coefficients of (±l) r W n (/3 — l)L™/ 2 J/[n/2j!, where [-J denotes the greatest integer less than or 
equal to its argument. 

Now, for any ip £ J?"(R), £±((3) := <p{Ep^±) is a holomorphic function of (3. To see this, we note that 

(11) £±(0) = V(±l)"e„AA» lP / = V(AA 2wC2 »±AA 2 » +1 c 2 » +1 ) lP , J , 

* — ' \n/Z ! * — ' n! 

?i=0 L ' J n=0 

where c„ are the Hermite coefficients of Lp. Since the c n are bounded by a polynomial in n, the series converges for 
all /3 £ C, giving holomorphy. Thus, if £±{/3) = for all /3 in an interval (as is the case here), then all of £±'s power 
series coefficients are zero. Applying this result to the two choices of sign, we obtain (since the Af n are never zero) 
c„ = V r) € No = 0, which thus proves the lemma. □ 

Remark. This result may be of wider applicability, particularly in harmonic analysis, due to the ubiquity of the 
Gaussian. We thus note that the proof of the lemma shows that a need merely belong to some subset of the right 
half-plane with an accumulation point to guarantee density. One could have also proved this result more abstractly 
(and without recourse to the Hermite expansion) by a slightly indirect application of the Stone-Weierstrass theorem, 
though the basic Hahn-Banach argument (contained in the Corollary from Conway we use) remains the sameQ 

5. Proof of TheoremQ] 

We first note that Eq. (O is clearly true for d = 1 (indeed, d E N) by the standard Poisson summation formula 
for lattices (applied to Z d ). To prove the result for d > 1, we shall first establish that it holds for the Gaussians from 
LemmaQ] and then show that the equality still holds in the limit in the Schwartz space topology. The control afforded 
by demanding convergence in the Schwartz space makes this quite straightforward. The primary result that needs to 
be shown is that two functions that are e-close in the Schwartz space topology have dimensionally continued Fourier 



Nota bene: Simon defines the h n without the factor of (— l) n (that here comes from our H n ). We have included the (—1)™ for notational 
simplicity (since we use the standard convention for the Hermite polynomials). This does not have any effect on Simon's Theorem 3, since it simply 
amounts to a sign change of the odd Hermite coefficients. 
6 This is 05.01.1 1.0001.01 on the Wolfram Functions Site BT1 . 
7 Personal communication from John Roe. 
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transforms that are Ce-close in a given Schwartz space seminorm (where the constant C depends on the seminorm 
under consideration, as well as d). 

To show that Eq. (O holds when / = Q a , we first consider the left-hand side and note that 



(12) 



-al 



0(e" 



Convergence is guaranteed because is analytic inside the unit disk. [To see that is analytic inside the unit disk, 
note that t?3 is analytic there, and, moreover, nonzero, so its dth power is analytic, as well. It is easiest to see that 1)3 
is nonzero inside the unit disk from its infinite product expansion, given in, e.g., Eq. (35) in Chap. 4 of Conway and 
Sloane ITOl .l Using LemmaQ] the right-hand side of Eq. (O becomes 



d/2 



1=0 

Now, the Jacobi imaginary transformation for ^3 [Eq. ©] implies that (n / a) d / 2 Q(e~ 7T l a ) = 0(e _Q ), so we have 
thus established the result for Q a . 

We shall now show that this equality continues to hold in the limit. The equality is clearly true for any finite 
linear combination of the Gaussians from Lemma [T] so we use Lemma [2] to approximate an arbitrary / € S fiE (M.) 
by a finite linear combination of these Gaussians, g. Specifically, we have j|/ — ,g|j Tl . m < e V n,m <E No, where 



l/ll 



su PzeR \x n f( m >(x)\ is the family of seminorms that gives the Schwartz space topology. (We denote the 



mth derivative of / by f^.) We wish to bound the difference between the two sides of Eq. (O by a constant times e. 
We have 



(13) 



OO 

E 

1=0 



00 



1=0 



< 



1=0 



1=0 



Ni(f-9)(Vl) 



where we used the fact that the dimensionally continued Poisson summation formula holds for g, along with the 
triangle inequality. We can bound the two sums on the right-hand side by constants times e using the assumption about 
the closeness of / to g in the Schwartz space topology and the fact that Ni grows at most polynomially with I. The 
latter fact also shows that the two sums on the left converge for / <G 

5.1. Bounds on the growth of 2Vj and on the right-hand side of Eq. ( fl3l i. 

Lemma 3. The Ni defined in Eq. d 1 Ob are polynomially bounded. Specifically, we have 

\Ni\ < 2 d (l + d/l) l (l + l/d) d <C d l d , 
where C' d > is some constant (and the second inequality only holds for I > 1). 

Proof. Recalling that is analytic inside the unit disk, we can apply Cauchy's integral formula to the contour Cr, a 
circle of radius R 6 (0, 1), centered at the origin (and oriented counterclockwise), to obtain 





f di{z) dz 


1 




2ni 




~ 2^ 


Jo 



R l e ilt 



-d9 



< 



R l {l-R) d ' 



Here we have used | ^3(9) | < 2/(1 — |g|) (for \q\ < 1, obtained using the geometric series). The right-hand side attains 
its minimum [for R G (0, 1)] at R = 1/(1 + d), so we have the desired result. To obtain the second inequality, we use 
the fact that (1 + 1 /r) r < e for r > 0. □ 

Remark. While this bound is easy to obtain and is all that is necessary for our purposes, it is by no means optimal. 
For instance, for integer d, we can apply the Hecke bound for modular forms (e.g., Theorem 4.5.3 in Rankin BTI ) to 
conclude that Ni = 0(l d ^ 2 ). (Other of Rankin's results — stated as Theorem 12 in Pache ll30l — show that the Hecke 



bound stated by Rankin is applicable to ^3 for d £ 



■) 



If we write h := f - g, then this bound implies that \Nih(Vl)\ < C d l d \h(y/l)\ < eC d /l 2 (for I > 1), where the 
second inequality follows from the fact that h is e-close to in the Schwartz space topology. [Explicitly, we have 
\x 2d+4 h(x)\ < e V x > 1 =*> l d \h(Vl)\ < e/l 2 V I G N. The first inequality comes from noticing that for any 7 > 0, 
we have |a; 7 /i(x)| < \x^~^ h(x)\ < e for x > 1, where ["•] denotes the smallest integer greater than or equal to its 
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argument.] We shall show that \p 2n h(p)\ < Kd£ V n £ N, p 6 R (where Kd is some rt-dependent constant), so we 
have \p 2d+i h(p)\ < Kdt V p <E R. We can thus apply the same argument to the second sum and hence bound both 
sums by constants times e (since Ya^i ^ 2 i s finite), showing that the dimensionally continued Poisson summation 
formula is true in the limit [since we will have shown that the right-hand side of Eq. (fT3l l is bounded by a constant 
times e]. 

5.2. Bound on \p 2n h(p)\. To prove the bound on \p 2n h(p)\, we first dimensionally continue some standard Fourier 
results. 

Lemma 4. If we define the d-dimensional Laplacianfor radial functions by 

(14) A d f(r) := f"(r) + —f'(r), 

r 

then, for d > 1, 

i) F p {r) := Fx(d/2; —n 2 p 2 r 2 ) satisfies A d F p = -A-K 2 p 2 F p , so 

ii) A«/(p) = (-ir(2Trp) 2n f(p)forfe S*(R). 

Proof. Part / follows from the fact that y a (r) := oFi(a; r) satisfies ry"(r) + ay' a (r) — y a (r)^ [Alternatively, it can 
be obtained by direct calculation using Eq. Q, justifying term-by-term differentiation using analyticity.] Part ii is 
then obtained by induction, applying Eq. (O to A' d 1-1 / and integrating by parts twice. [The boundary terms at infinity 
vanish because / E 5f (R); those at vanish because d > 1 (or cancel amongst themselves).] □ 

We can thus write \p 2n h(p)\ = (27r)~ 2n |A^ft,(p)|. Then, since we shall show below that \r k A2h(r)\ < T>e, where 
D is some (n- and d-dependent constant), we obtain [using Eq. © and the fact that oFi(a; r) is a bounded function of 
r, as was seen in the proof of LemmaQ") 

/>oo 

\p 2n h{p)\ <C I lAS^rJIr^dr 

f 1 poo 







<C 

<cv 



\A n d h(r)\dr+ lA^WIr^^r 
o J\ 



s dr 



where C > is some (n- and d-dependent) constant and we used \r k A^h{r)\ < Ve with k = and k = s. We can 
choose s > d, so the integral in the final term is finite, thus giving the desired result. 

5.3. Bound on \r k A^h(r)\. To see that \r k A^h{r)\ is bounded by some (n- and (^-dependent) constant (called V 
above), we first note that we can use induction to write 



(15) A3fc(r) = $> 



2n fcC*)( 



r 



1 r 2n—j 
3=1 

for some (n- and d-dependent) constants aj (and an arbitrary differentiable/i). Thus, for |r| > 1, we have \r k A^h(r)\ < 
e Y^j=i a j- For \r\ < 1, matters are considerably more subtle, and we have to rely on the fact that h is even to see 
that Ajh remains bounded at the origin. The argument goes as follows: We write h = V + 7Z, where V is h's 
(2n)th-degree Maclaurin polynomial (necessarily even, since h is) and 7Z is the associated remainder. We then have 
\r k A^h{r)\ < \r k A^V (r)\ + \r k A^ l 72.(r)|. Since AJJ maps even polynomials to even polynomials [as can be seen 
from Eq. (fTBIll, \r k A^V (r)\ is bounded by a (fc-, n-, and d-dependent) constant times e for |r| < 1. [Since the co- 
efficients of V are given by derivatives of h, they are bounded by constants times e, by hypothesis.] To deal with 
|r fe A27?.(r)|, we first need to establish an identity for derivatives of 7Z, viz., (for j < 2ri) 



[) (2n + l-j)! 



r 2n+l-j 



8 This differential equation for -Fi is 07.17.13.0003.01 on the Wolfram Functions Site ED. 
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for some £j € (0, r). One obtains this by comparing the jth derivative of h = V + 1Z with the (2n — j)th order 
Maclaurin expansion (with Lagrange remainder) of hp). The polynomial pieces are the same, while the remainder 
pieces give the two sides of the equality. Combining this identity with Eq. 1151 . we obtain 

2n 

i=i 

where the bj are (n- and <i-dependent) constants. This shows that \r k A^lZ(r)\ is bounded by an (n- and d-dependent) 
constant times e for \r\ < 1, so \r k A~2h(r)\ is, as well, proving the desired result, and hence the theorem. □ 

Remark. The restrictions on / and d in the statement of the theorem are surely not optimal: Bochner's proof of the 
basic summation formula holds for d € M, d > 0, and / of only finite differentiability (see Theorem 9 in 0), and 
there is numerical evidence that the given result holds for d 6 C, Red > and significantly less smooth / [e.g., 
f(r) = e~' r l , which is not covered by Bochner's results]. (The evidence also extends to the generalization given in 
Theorem[2]and is provided by a MATHEMATICA notebook, available onlinefl) While one could use a slightly larger 
function space than ,5^ E (R) without any change to the proof — the proof does not need control over ||/ — <7|| ra ,m for 
all n and m — we did not investigate this in any detail: The resulting function space would still require a fair amount 
of differentiability (while we have numerical evidence that the formula remains true for at least some functions with a 
cusp at the origin), and faster decay than the standard Poisson summation formula (see, e.g., Corollary 2.6 in Chap. VII 
of Stein and Weiss ll37l ). Moreover, the closure of the family of Gaussians in this less restrictive topology would almost 
surely be more recondite than <9' E (M.). 

6. Generalization of Theorem[T] 

Since there are other families of lattices with dimensionally continued theta series besides Z d (e.g., the root lattice 
D d mentioned in Sec. 12.1b . it is reasonable to expect that Theorem [TJ can be generalized by replacing with some 
more general function T, which we shall term a generalized theta series. It is not clear how to construct the most 
general such T0 However, the template provided by the theta series of other standard lattices (e.g., the ones given 
in Chap. 4 of Conway and Sloane 1101 ) allows us to construct a reasonably general family of generalized theta series 
out of finite linear combinations of products of the three Jacobi theta functions ($2, $3, and $4) given in Eq. (f3]l. This 
family will contain all of the theta series of lattices and shifted lattices given in Chap. 4 of Conway and Sloane iflOl . 
except for the general form of the theta series of the root lattice A d and its translates. In fact, theorems in Conway and 
Sloane (Theorems 7, 15, and 17 in Chap. 7 and Theorem 5 in Chap. 8) show that the theta series of large classes of 
lattices can be written in such a form. However, the family of Ts is considerably more general, since one only requires 
A m , Pm, &m G No to reproduce the theta series in Conway and Sloane, while here they can be arbitrary nonnegative 
real numbers. The Ts also contain several families of modular forms, as shown by Rankin OTl . and concisely stated 
in Theorem 12 of Pache BUI . 

Explicitly, we make the following 

Definition. A generalized theta series is a finite linear combinations of T,js of the form 

M 

(16) T d (q) := ft 4-(q^)^(q^)^(q Um ), 

with \ m ,p m ,cr m > 0, J2m=i(^n + Pm + & m ) = d, and s m , t m , u m G Q+. 
We thus have 

_ £ *-tf'-*r<,v-)«i-(.v-) t 

x \ / A m ,p m o-,„ 

m—1 V °m b m u m 



^The notebook is available at |http : / / gravity ■ psu ■ edu/ -nathan jm/Dim_cont_PSF_test ■ nb| 

'"But note that Ryavec characterizes all admissible Ts (under certain assumptions) for d = 1 in 1331 . We also call attention to the work of 
Cordoba I11II12I . who shows that in integer dimensions, large classes of generalized Poisson summation formulae arise from the standard Poisson 
summation formula applied to the finite disjoint union of (integer dimensional) lattices. (Note that Lagarias makes a slight correction to the statement 
of Theorem 2 of 1111 in Theorem 3.7 of 1231 .) 
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which we compute using 

T*(z) := (z/z) d / 2 T(-lA). 

Nota bene: The definition of T* is just the dimensionally continued Jacobi transformation [Eq. ©J of T with the 
factor of vdetT omitted. (We leave off this factor, since it would just cancel against the one present in the standard 
Poisson summation formula [cf. Eq. ©].) 

We now need to show that the generalized theta series have the appropriate properties to allow us to copy the proof 
of TheoremQ] almost verbatim. Specifically, we need the following 

Proposition. IfT denotes any of the T^s or T* d s defined above, we can write 

oo 

1=0 

where 

1. A l+1 > At, A > 0. 

3. There exists LeN and C,n > such that \ Ni \ < CAf for all I > L. 

4. The series converges inside the unit disk. 

Proof. First note that we can apply the same arguments to T* d as to T^, so we can restrict our attention to the former, 
without loss of generality. Now, we then have 

M 

A t = {l + A)/V, A:=J2 VX rnS m /4: 

m— 1 

where V is the least common denominator of s m , t m , and u m (for all to). [We have the additive constant A due to the 
overall factor of q 1 / 4 in $2(3)-] Thus the first required and second required properties (positivity and monotonicity of 
the A{ and convergence of the series whose terms are A~[ 2 ) are obviously true. 

For the third property (polynomial boundedness of the Ni), we use the same Cauchy's integral formula argument 
used in Sec. 15. II (The analyticity established below shows that Cauchy's theorem is still applicable here.) Here Ni is 
given by the Zth term in the Maclaurin expansion of Td{q v )/q A , so we have 



\Ni\ 



1 f ft ^(^)^(^ tm )^'"(^"") , 

11 ~,V\ m s m /4.,l+l aZ 



2-Kl c„ 11 z V\ mSm /A z l+l 
H m— 1 

M 



<2 d \ 



< 



i i? ; (i-i? v/Sm ) Am (i-^ 1/t " l ) pm (l-^ VUm ) ,Tm " R l {i-R) d ' 

where Cr is the same contour used previously. We have used the geometric series to obtain the bound \r(q)\ < 
2/(1 — \q\), where r(q) is any of $ 2 {q)/q 1,A , tiz{q), or$ 4 (q). Additionally, we have used the fact that k > 1, where 
k is any of Vs m , Vt m , or Vu m , so |1 — _R K | > 1 — R, since R e (0, 1). We also recalled that A m , p m , a m > and 
Em=i( A m + p m + cr m ) = d. Since there is an R 6 (0, 1) such that 2 d /[R l {\ - R) d ] < C d l d (for I > 1), as was 
shown in Sec. 15. II we are done. 

The fourth property (convergence of the g-series in the unit disk) follows from the analyticity and lack of zeros 
of the theta functions inside the unit disk, as in the argument given below Eq. (fT2l . [Note that here we consider 
#2(<?)/Y /4 , not § 2 {q) itself.] Specifically, T{q v )/q A is an analytic function of q inside the unit disk; the lack of 
zeros can be seen from the infinite product representations of $2 an d ^4 given, e.g., in Eqs. (34) and (36) in Chap. 4 
of Conway and Sloane ifTOl . □ 

The dimensionally continued summation formula for generalized theta series thus takes the form of the following 

Theorem 2. Let T be a generalized theta series (as defined above), with power series coefficients and powers iVj and 
A h i.e., 

00 

T(q)=J2Niq A '. 

1=0 
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Let and A* be the corresponding quantities for T 's dual, T* [computed using Eq. d 1 7b ]. Then, for any f £ , 
we have the summation formula 

oo oo 

£jv I /(^)=x;jv?/(v^) ) 

(=0 (=0 

where we compute f using Eq. §6$ (with the dimension parameter d associated with T). 

Remark. In general, all one requires of the Ts used in this summation formula is that they and their Jacobi transfor- 
mations have sufficiently well-behaved power series. (This is satisfied by all functions of weight, in the terminology 
of Pache l30l — see Pache's Definition 11.) However, it is unclear whether any such functions exist besides the gener- 
alized theta series we have defined above, so we have not given the theorem in a more general form. 

Proof. The proof is almost the same as that for Theorem Q] (replacing by T, and noting that we can no longer 
appeal to the standard Poisson summation formula for d = 1, so we simply exclude that case). Most of the work 
has been done in the proof of the Proposition; the only new part is checking that J^hLo \Nih(\/A~i)\ —> as e — > 
if h is e-close to in the Schwartz space topology [and similarly for Yli^Lo l-^z*' l (\/^f)|]- To do this, we simply 
note that we have \Nih(y/A~i)\ < CA?\h(y/A~i)\, by polynomial boundedness of the N t , and that x 2n+i \h(x)\ < e 
V x E R => 1 7^ ( | < [cf. the discussion at the end of Sec. 15.11 , from which the desired result follows 

immediately. (The same argument holds for the starred quantities, since they have the same properties as the unstarred 
quantities.) □ 

Remark. This theorem can likely be interpreted as a trace formula for the dimensionally continued, spherically 
symmetric Laplacian [Eq. ([Pflil. since the kernel of the dimensionally continued Fourier transform is an eigenfunction 
of this operator (see Lemma|4|i. See, e.g., Sec. 1.3 (particularly Theorem 1.3) of Uribe l39l for a presentation of the 
standard Poisson summation formula for an integer dimension lattice as a trace formula for the Laplacian. 

Remark. This result shows that one can apply this extended Poisson summation formula to lattice-like objects whose 
theta series have coefficients of both signs, so they do not exist as a lattice, even though d € N: For a trivial example, 
consider d = 2 and Q(q) = #4 (q) = 1 — 4<j + Aq 2 + ■ ■ ■ . Of course, this is in some sense too trivial, since one can write 
#4 = 20 D 2 — 0^2, and then apply the standard Poisson summation formula to each of those lattices to establish the 
result in this case (cf. the discussion in Cordoba lfl2l ). However, in more complicated higher-dimensional cases, it will 
likely not be clear how to construct the lattice(s) associated with the theta series (if they indeed exist). Indeed, Jenkins 
and Rouse fl9l have very recently shown that for weights higher than 81 632, all the modular forms of a certain type 
have coefficients of both signs. (This general property was first shown in less specific, much earlier work by Mallows, 
Odlyzko, and Sloane 11251 .) 

7. Outlook 

While Theorem [2] encompasses quite a large family of summation formulae, there still remains wide latitude for 
further generalizations (even excluding the various possibilities for weakening certain of the hypotheses mentioned 
after TheoremQ]). The most sweeping generalization would likely be to replace the dimensionally continued Fourier 
transform with some more general family of integral transforms, with the possibility of a subsequent enlargement of 
the transformation properties required of the generalized theta series. Here one could follow the work of Kubota [22] 
and Unterberger ll38l in integer dimensions. But even if one retains the dimensionally continued Fourier transform, 
one can still likely obtain summation formulae from more general classes of generalized theta series than we have 
considered. In particular, it would be interesting to obtain a dimensionally continued version of the quasicrystal 
summation formula given as Theorem 2.9 in Lagarias (23). Here the calculations of the central shelling for certain 
quasicrystals in, e.g., [3] |27l [28J |291 could be relevant. Additionally, since the coefficients of standard theta series 
give the representation numbers for lattices, it is possible that our results could be applicable to generalizations of 
representation number problems: See, e.g., |[T6l for a review of standard results on representation numbers. 
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